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Jamming V: Jamming in two-dimensional packings 
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Abstract. - We investigate the existence of random close and random loose packing limits in two- 
dimensional packings of monodisperse hard disks. A statistical mechanics approach — based on 
several approximations to predict the probability distribution of volumes — suggests the existence 
of the limiting densities of the jammed packings according to their coordination number and 
compactivity. This result has implications for the understanding of disordered states in the disk 
packing problem as well as the existence of a putative glass transition in two dimensional systems. 
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Introduction. — The concept of jamming is a com- 
mon feature of out of equilibrium systems experiencing a 
dynamical arrest ranging from emulsions, colloids, glasses 
and spin glasses, as well as granular materials [1]. For 
granular matter, it is argued that a statistical mechanical 
description can be used, with volume replacing energy as 
the conservative quantity [2]. In this framework, a meso- 
scopic model has been presented [3] , allowing the develop- 
ment of a thermodynamics for jamming in any dimension. 
Here we develop this theoretical approach to investigate 
■ the existence of disordered packings in two-dimensional 
systems composed of equal-sized hard disks [4] . The exis- 
tence of amorphous packings in 2d is a problem of debate 
in the literature: two dimensional systems are found to 
crystallize very easily since disordered packings of disks 
are particularly unstable [5] . 

In two dimensional Euclidean space, the hexagonal 
packing arrangement of circles (honeycomb circle packing) 
' has the highest density of all possible plane packings (or- 
dered or disordered) with a volume fraction <j>hex = — 
0.9069 and each disk surrounded by 6 disks. Regarding 
amorphous packings, experiments find a maximum density 
of random close packing (RCP) of monodisperse spheres 
at rcp « 0.82 [4] while the lower limit (random loose 
packing, RLP) has been little investigated and its exis- 
tence has not been treated so far. The tendency of 2d 
packings to easily crystallize has led to consider bidisperse 
systems which pack at a higher RCP volume fractions of 
rcp «O.84 [6,7]. 

In parallel to studies in the field of jamming — which 
consider the packing problem as a jamming transition ap- 
proached from the solid phase [6-8] — other studies at- 



tempt to characterize jamming approaching the transition 
from the liquid phase [9, 10]. Here, amorphous jammed 
packings are seen as infinite pressure glassy states [9]. 
Therefore, the existence of disordered jammed structures 
(of frictionless particles) is related to the existence of a 
glass transition in 2d [9] , a problem that has been debated 
recently [11]. 

Here we treat the disordered disk packing problem with 
the statistical mechanics of granular jammed matter [2]. 
The formal analogy with classical statistical mechanics is 
the following: the microcanonical ensemble, defined by all 
microstates with fixed energy, is replaced by the ensemble 
of all jammed microstates with fixed volume. Hence, the 
appropriate function for the description of the system is 
no longer the Hamiltonian, but the volume function, Wi, 
giving the volume available to each particle unit such that 
the total system volume is V — J2i W« P]- 

The aim of the present work is to develop the model 
presented in [3] for the calculation of the volume function 
in the case of 2d packings. The validity of the hypothe- 
sis employed in [3] are discussed, and they are modified 
according to the properties of 2d packings. We use our 
results to study the nature of the RLP and RCP limit 
in 2d through an elementary construction of a statistical 
mechanics that allows the study of the existence of a max- 
imum and minimum attainable density of disordered circle 
packings. We find that amorphous packings can pack be- 
tween the density limits of ^77.5% and ^80.6% defining 
the RLP and RCP respectively, according to system co- 
ordination number and friction, opening such predictions 
to experimental and computational investigation. While 
these values should be considered as bounds to the real 
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values due to the approximations used in the theory, they 
serve to suggest the existence of both limits in two dimen- 
sional packings of monodisperse disks. 

Volume function. — The volume function is the key 
of the system's statistics: its flat average over all jammed 
configurations determines the total volume. The most nat- 
ural way of dividing the system is called the Voronoi dia- 
gram, which can be seen in Fig. la. Each grain's region 
is the part of the space closer to this grain than to any 
other, so that the volume is clearly additively partitioned. 
The major drawback of this construction is that, so far, 
there was no analytical formula for the Voronoi volume of 
each cell, such that attempts have been made to use other 
constructions [12]. 

The Voronoi volume of a particle i can be written as: 
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where j\j is the vector from the position of particle i to 
j, the average is over all the directions s forming an angle 
0ij with fij as in Fig. la. This formula has a simple inter- 
pretation depicted in Fig. la. For consistency of notation 
with previous work, we will use the words "volume" and 
"surface" in 2d, although they correspond to "surface" 
and "length" respectively. 

The volume function defined in terms of the particle 
coordinates is of no use, since it does not permit the cal- 
culation of the partition function. To solve this problem, 
we calculate an average free volume function based on the 
environment of the particle, referring to a coarse-graining 
over a certain mesoscopic length scale. We assume a prob- 
ability distribution for the positions of the nearest neigh- 
bors as well as the other particles. After averaging over 
the probability distribution we obtain an average meso- 
scopic free volume function representing quasiparticles in 
the partition function. Considering isotropic amorphous 
packings allows for removal of the oricntational averaging. 

Probability distribution of volumes. Using the 
notation of Fig. lb, we see that the microscopic vol- 
ume function is entirely defined by the parameter c = 
min[r/ cos#]. The calculation of the average free volume 
function, w, requires knowledge of the probability distri- 
bution of this parameter. That is: 



w ee (wr)/v g 



/■OC 

-1 = -/ (c 2 - 

Jc=l 



1) 



dc 



dc, (2) 



where P>(c) represents the inverse cumulative distribu- 
tion, i.e. the probability that all balls verify r^/ cosOij > 
c, which is calculated under the following hypothesis: 

(1) The cumulative distribution P>(c) is made of two 
contributions, one of the "contact" balls (i.e. touching 
the considered grain), called P> , and one of the other 
(background) balls, P> . These probabilities can be un- 
derstood as the probabilities of a particle in contact (resp. 
background) for being situated outside the grey zone on 
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Fig. 1: (a) The limit of the Voronoi cell of particle i in the 
direction a is rij/2 cos Oij . Then the Voronoi volume is propor- 
tional to the integration of (rij/2 cos Oij) over s as in Eq. (1). 
(b) The particle contributing to the Voronoi volume along s is 
located at (r,ff). The dark gray region is the considered grain 
(r < R), and in white the excluded zone for the center of any 
other grain (r < 2R). For a given c = r/cos#, the light grey 
area is the region of the plane (r',6 1 ) where r' / cos 9' < c. 



Fig. lb, and therefore not contributing to the Voronoi 
volume c 2 . 

(2) The probability distributions are those of a large 
number of particles at a given density, and of negligible 
size giving rise to Boltzmann-like distributions: 



P> B (c) 



exp(- P V*(c)) 



and 



P°(c)=exp(-p s Sr(c). 



Here V* and S* represent a free volume and surface re- 
spectively towards c, i.e. V*(c) = J 6(c — j^~)dr and 
S*(c) = § 0(c — j7g)ds where the integrals cover respec- 
tively the space and the unity sphere. The densities, p(w) 
and ps(z), are mean free- volume and free-surface densi- 
ties, respectively. 

In 2d, the volume of the grain (with 2R = 1) is V g = 
7r/4. The free volume density (inverse of the free volume 
per particle) is p(w) = N/(NV g (t)- 1 - NV g ) = l/(V g w). 
Then, 



V*(c) = 



- 1 



arccos(l/c) 



1 - 



,2 ' 



The main assumption here is that the packing structure 
is uniform, thus the pair distribution function is assumed 
to be a delta function at contact plus a constant for larger 
distances. This assumption is an oversimplification, and 
more realistic background could be considered, such as 
peaks in the pair distribution at the next nearest neighbor 
sites. 
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For the surface contribution, wc have: 



/>arccos(l/ c) 

S*(c) =2 d9 
Jo 



2 arccos(l/c) 



The surface density ps = 1/ (S) is the inverse of the av- 
erage surface left free by z contact balls (see Fig. 2a). As a 
rough approximation, one can assume it is proportional to 
z, but because of the size of one ball, there is an "excluded- 
surface" effect, so that the exact value is determined by 
numerical simulations. It consists in setting sequentially 
and randomly z non-overlapping circles of radius 1 at the 
surface of the unity circle (Fig. 2a). The closest ball to 
the considered direction s defines the free angle. The free 
surface is then twice this angle. Its average value is the 
mean free-surface (S). 
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Fig. 2: (a) Free surface, in an example of z — 4 contact balls, 
denned by the angle of the closest grain to s, with (S) = 26* . 
(b) Simulation results for 1 < z < 5, with a second-order 
polynomial fitting, and linear fitting for 3 < z < 4 : ps ~ 
(z-Q.5)/n 

Results are shown in Fig. 2b. Important deviations 
from the linearity in z are not surprising, since each con- 
tact ball occupies an important surface (z max — 6 in 2d), 
and strong finite-surface effects are expected. For z = 5, 
in around 41% of the trials, the fifth ball cannot be set be- 
cause there is not enough space, and we take into account 
only the 59% remaining trials. 

For 3 < z < 4, we will use the linear dependency 
Ps(z) = 2 ~ 7r °' 5 as fitted in Fig. 2. Obviously a fitting for 
1 < z < 5 would be of higher order, but in our range, the 
error is insignificant compared to other approximations of 
the model. 

(3) The cumulative distributions are not independent. 
The assumption that the surface and volume terms do not 
overlap seems to be an abusive approximation in 2d. This 
is not the case in 3d as shown in [3]. Indeed, for higher 



dimensions the probabilities are expected to become inde- 
pendent, but in the case of 2d a new solution has to be 
worked out which considers the correlations between the 
contact and background term. 

To illustrate this point, Fig. 3a shows the considered 
grain, the free volume V*(c) and the circle occupied by 
the closest "contact-grain" (i.e. the excluded zone for the 
center of any other grain because of its presence) for a 
value of c — 1.2. In fact, the values of c < 1.2 are con- 
tributing for 94% of the distribution if we neglect the over- 
lap of contact and background grains. As we see, the free 
volume is mostly covered up by this contact-grain, and the 
non-overlapping hypothesis used in [3] appears obviously 
wrong. This statement is confirmed by the calculation 
of the RCP density: with the non-overlapping hypothe- 
sis, the calculation provides a value of 4>rcp = 0.89, to 
be compared with the reported value of 0.82. The nearer 
grains are exceedingly taken into account. 
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Fig. 3: (a) For c = 1.2, the star is the center of the clos- 
est contact-grain, occupying the circular region printed with a 
pattern. The free volume is printed in grey. It is almost com- 
pletely overlapped by the surface contribution, (b) The closest 
contact-ball depends on cc, the free volume on cb, and AV 
on both. Here it is the intersection between the region in light 
grey and the patterned region. 

Therefore, the volume term is modified by substitut- 
ing the free volume V* by V* — AV* which represents 
the free volume minus the part occupied by the closest 
surface grain. The meaning of this change is that the 
contributions are no longer independent, and depend on 
two parameters cb and cc- The distribution P>(c) is the 
probability that both cc and cb be higher than c. Figure 
3b shows the overlap of the contact grain parameterized 
by cc and the background volume parameterized by cb, 
defining AV* (cb, cc). The analytical formula of AV* is 
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determined by geometrical calculations. 
The probability density is 

P(c) = -dPy/dc = -P°(c) ■ dP^/dc - P*(c) ■ dpC/dc. 

The meaning of the latter equality is that, to realize c, we 
must have either cc or cb equal to c, and the other higher. 
The background probability depends on cq'. 




P*(c B \c c ) = exp[-p(l/*(c B ) - Ar(c B ,cc))], 



(a) 



and 



P B (cb\cc) = -^P>(cb\c c )- 

CICb 



Then, 



P(c) 



P c (c c )P B (c\c c )dcc 



P C (c)P*(c\c) 




d_ 

dc 



P (c c )P^(c\c c )dc c = -^P>(c). 



(b) 
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From (2), we integrate by parts using the latter equal- 
ity. The boundary term [(c 2 — l)P>(c)] vanishes, since 
the limits of integration correspond to c = 1 and c — > oo, 
with P(c — > oo) = 0. We obtain for the average volume 
function from Eq. (2): 



w 



POO PC 

> C 

Jc=l Jc 



dpv 

dc c 



with P^(c c ) 



exp[-pM(V*{c)-AV*(c, c c ))]dc c dc^ 

(3) 
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(c) 



, - — exp[— ps(z)S*(cc)]- Fig. 4: (a) w(z) curve, with a linear fitting: w(z) = 0.437 — 

0.0492. (b) Volume fraction according to the second Eq. (4). 

Free volume function. - Equation (3) is a self- (c) Prediction of the model The thick curveg are the limit of 

consistent equation to obtain w(z), which cannot be solved the diagram at x = and X -> oo. We show several curves 

exactly, therefore a numerical integration of (3) is neces- G f constant compactivity, X. The curves are plotted for (from 

sary to obtain w vs z. For various values of z, we integrate right to left) : X = 5.10~ 3 , X = 10~ 2 and X = 10 _1 . The 

Eq. (3) numerically, and we then calculate a fitting of the horizontal lines are both constant Z given by an arbitrary no 

results (Fig. 4a). We obtain the free volume function and (dotted) and — > oo (thick inferior limit of the diagram). 

the local density cj)' 1 = w + 1 (Fig. 4b): 



w(z) = 0.437- 0.049z, <j>i(z) 



1 



1.437-0.049^ 



(4) 



Generally speaking, we would expect w to be roughly 
proportional to 1/z, with w — > when z — > oo. However, 
the statement z — > oo has little meaning when we plot a 
figure for 3 < z < 4 and the "infinite" (maximal) value of 
z is 6. 

Statistical mechanics. — Equation (4) plays the role 
of a " Hamiltonian" of the system. Each jammed configu- 
ration corresponds to some " volume level" in analogy with 
energy levels in Hamiltonian systems. From the formal 
analogy with classical formulas, the canonical partition 
function is [2]: 



Z(Z,X) = J g{w)e- w/x dw, 



(5) 



where X is the (reduced) compactivity (normalized by the 
volume of the spheres) and g(w) is density of states for a 



given volume w. We remind that the compactivity is the 
equivalent of temperature in the Edwards statistics, and 
it is a measure of the system's looseness. 

Since the volume w is now directly related to coordina- 
tion number z through Eq. (4), we can compute g{w) by 
replacing variable, g(w) = J P(w\z)g(z)dz, where P(w\z) 
is the conditional probability and g(z) is the density of 
states for given z [3]. 

At this point a distinction has to be emphasized: we 
refer to z as the geometrical coordination number since it 
is purely defined by the particle positions. On the other 
hand, there is the mechanical coordination number, Z, 
defined by those geometrical contacts that carry a non- 
zero force. Z is then defined by the mechanical constraint 
leading to the isostatic condition. That is, Z = 2d = 4 for 
frictionlcss disks while Z = d + 1 = 3 for infinitely rough 
disks. For a system with a finite interparticle friction fi, 
Z(p) interpolates between both limits [3,7]. 

Obviously, z must be larger than Z for the mechani- 
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cal condition to be satisfied. The different volume lev- 
els of a packing can be understood in the following way: 
the friction coefficient sets a mechanical constraint on Z, 
but the system can explore all geometrical levels z > Z. 
Additionally z is bounded by the maximal coordination 
number for a random packing, which is 2d = 4, since 
there are z/2 constraints on the d particle coordinates. 
In relation with the discussion of isostaticity, it is believed 
that above this value, the system is partially crystallized. 
Therefore we obtain: g(w) = J z P(w\z)g(z)dz. Since w(z) 
from Eq. (4) is a coarse-grained free volume and inde- 
pendent of the microscopic positions of particles, we have 
P(w\z) = 8(w — w(z)). 

The expression of the density of states is obtained 
by considering that the states are collectively jammed. 
Therefore, the space of configuration is discrete since we 
cannot continuously obtain one configuration from an- 
other. Assuming a typical distance between configurations 
as h z , we obtain g(z) cx (h z ) z , the exponent z arising since 
there arc z position constraints per particle in the jammed 
state compared to the free (gas) state. Such a formula is 
analogous to the factor h~ d for the density of states in tra- 
ditional statistical mechanics, where h is the Planck con- 
stant, which arises because of the uncertainty principle, 
i. e. because of the discreteness of the elementary volume 
of phase space. 

Substituting into Eq. (5), we get : 



,4 

Z(Z,X) = / (h z ) z e~ w( - z ^ x dz 
Jz 



(6) 



To establish the maximum and minimum densities, we 
consider the limits of zero and infinite compactivity, re- 
spectively. The ground state of jammed matter, is anal- 
ogous with the limit T — > 0. The only accessible state is 
z = 4, corresponding to the random close packing. For this 
state, from Eq. (4) we get a fixed value of the volume frac- 
tion for any coordination number Z G [3,4], /j 6 [0, oo]: 



y rcp 



(Z) 



1 



1.437- 0.049 x 4 



0.806, 



(7) 



The lower density appears for X — > oo when the Boltz- 
mann factor is unity in Eq. (6), and we obtain the densi- 
ties of RLP (assuming d z « 1): 



Z{Z, oo\ 



437 - 0.049z 



(h z ) z dz 



1 



Ze [3,4]. 



(8) 



1.437 -0.049Z' 

This leads to the diagram in the plane (</>, Z) plotted in 
Fig. 4c defining the possible jammed configurations. On 
the right part of the vertical line defined by Eq. (7), no 
disordered packing can exist. To the left of Eq. (8), the 
packings arc not mechanically stable. 

Between these two lines, we plot the lines of constant 
finite compactivity. For a finite value of the compactivity, 



the equation <j>{Z) is calculated by numerical integration. 
In the figure, three curves are plotted, respectively X = 
5.10~ 3 , X = IQ- 2 and X = 10" 1 for h z = 0.01. The 
compactivity increases from the right (X = 0) to the left 
(X — > oo). The limit /j, — > oo, Z — > 3 defines the lowest 
RLP density value which is predicted: 



1 



0.049 x 3 



0.775 (9) 



The value of r i p depends on the mechanical coordina- 
tion number, contrarily to the value of <f) lcp . The shape 
of the diagram is similar in 3d [3], and this is in agree- 
ment with the wide range of reported values for RLP, in 
contrast with RCP [4]. On a horizontal line given by a 
system with fixed Z, packings of different volume frac- 
tions can be achieved by applying different quench rates 
or compression speeds during the preparation protocol. 
Slow compressions achieve loose packings (and high com- 
pactivities). The obtained predictions for the density of 
RCP are close to the experimental values while we predict 
the existence of a RLP density. 

Summary. — In summary, we have used a model of 
volume fluctuations to develop a statistical mechanics of 
granular matter in 2d. From a quantitative point of view, 
we have seen how it lies on several approximations, that 
can appear too rough. The main difference with the 3d 
case is the need of taking into account properly the corre- 
lations in the probability distribution of volume through 
the consideration of point (3) above. Indeed, if we take 
Py and P> to be independent as considered in [3] we 
find r ip = 0.84 and rcp = 0.89, both values above the 
experimental value of RCP. The results, although not al- 
lowing exact predictions, are situated in the right order of 
magnitude for the limiting volume fractions. Due to the 
several approximations of the theory, the resulting limit- 
ing densities have to be considered as bounds to the real 
values. Improvements can be achieved by taking into ac- 
count the size and shape of the disks, as well as exact enu- 
meration to calculate P>{c), which can be done at least 
to a prescribed coordination shell of particles, in a brute 
force analysis analogous to the Hales proof of the Kepler 
conjecture, currently being undertaken. Altogether, the 
present framework seems to be successful in describing at 
least qualitatively the general features of jammed gran- 
ular matter in 2d providing evidence of the existence of 
RCP and RLP and their density value. These results sug- 
gest that a putative ideal glass transition may also exist 
in frictionless hard disk as discussed in [9] . 
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